On Upper Bounds for the Tail Distribution of 
Geometric Sums of Subexponential Random 

Variables 

Andrew Richards 



Department of Actuarial Mathematics and Statistics and the Maxwell Institute for 
Mathematical Sciences, Heriot-Watt University, Edinburgh 
E-mail awr2@hw.ac.uk 

The approach used by Kalashnikov and Tsitsiashvih for constructing up- 
per bounds for the tail distribution of a geometric sum with subexponential 
summands is reconsidered. By expressing the problem in a more probabilistic 
light, several improvements and one correction are made, which enables the con- 
structed bound to be significantly tighter. Several examples are given, showing 
how to implement the theoretical result. 

Keywords: Geometric Sum, subexponential distribution, upper bounds, GI/GI/1 

1 Introduction 

This paper was motivated by a series of papers by Kalashnikov and Tsitsiashvili [HI 
[19] . In these papers they proposed a method for constructing bounds for the relative 
accuracy of the asymptotic approximation to the tail distribution of a geometric sum. 
In trying to understand their method we found an error in their derivation which 
led to an incorrect formula (Lemma 1 in [18] and Lemma 6 in [19]). However, their 
general methodology was sound. We used a more probabilistic approach to deriving 
equivalent corrected results. This has simplified the picture, and the more natural 
probabilistic interpretation has enabled us to achieve much tighter bounds over a 
greater domain of application. 

Let Xi, X2, . . . be a sequence of non-negative independent and identically distributed 
random variables with unbounded support on the positive half-line. Let v be an 
independent counting random variable with geometric distribution, 

P(z/ = k)= p(l - pf'^, A; > 1, < p < 1. 

Then By = Yli=i -^i ^ Geometric Sum. We define 5*0 = 0. 

Many situations can be modelled as a geometric sum. Applications include risk pro- 
cesses, ruin theory, queueing theory and reliability models. The following list of 
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references contains useful surveys: Abate, Choudhury and Whitt p[], Asmussen [2], 
Cocozza-Thivent and Kalashnikov Feller [10], Gnedenko and Korolev [13], Gne- 
denko and Ushakov [T^ , Grandell [15] , Kalashnikov [lEl [17] . 

In queueing theory, there are many popular models where geometric sums arise natu- 
rally. The first classical example is for a GI/GI/1 queue with stationary waiting time. 
It is well known (see [TD]) that the waiting time coincides in distribution with the 
supremum of a random walk with negative drift, M = sup„ Yl^=i 65 where the are 
functions of the interarrival times and service times, and have common distribution 
function G{x). The supremum of a random walk can then be modelled as a geometric 
sum P(M > x) = {1- p)P{S^ > x), where = ^^"^^ ^pi and p = P(M = 0). The 
random variables ipi appearing in the geometric sum have the distribution of the first 
overshoot of the random walk over level zero, conditional on such an overshoot oc- 
curring. The distribution of the ipi does not coincide with G. However we can obtain 
bounds on the distribution of the ipi, (see, for example, Borovkov [5]). In [5j, Chapter 
4, Section 22, Theorem 10 we find 

ciG^(x) < P(V^i > x) < ciG^(x) + C2G(x), (1) 

where G{x) = 1 — G{x) and G^{x) = G{y)dy, and Ci and C2 may be estimated 

under further assumptions. We shall give an example of this in section 4. 
The tail distribution of the geometric sum is the object of interest. In the light-tailed 
case the situation is straightforward and an upper bound for the tail distribution is 
well known, the Cramer Lundberg upper bound. However, when the summands have 
a subexponential distribution it is much more difficult to find an appropriate upper 
bound for the distribution. Analytical techniques, simulation and recursive techniques 
have been used to estimate the tail. The tail is usually approximated by its asymptotic 
form (see, for example, Embrechts and Veraverbeke [9]), but this approximation can 
be very poor, and therefore it is of use to have an upper bound for the relative 
accuracy of the asymptotic approximation. A complete list of references for works 
on the bounds of the tail of the geometric sum distribution would be too long, but 
among important works are: Asmussen, Binswanger and Hojgaard [3j, Dickson and 
Waters [8], Kalashnikov [IH HZ] , Lin [21] , Willmot [22], Willmot and Lin [23]. 

The rest of this paper is structured as follows. In section 2 we lay the ground work and 
derive an auxiliary result concerning subexponential distributions. In section 3 we ob- 
tain the required upper bound for the relative accuracy of the asymptotic expression. 
In section 4 we apply the result to a number of subexponential distributions. 

2 Some Preliminary Results 

We consider distributions F{x) with unbounded support on the positive half-line. 
We will write the tail distribution as F{x) = 1 — F{x), and F*"(x) as the n-fold 
convolution. 

We recall the following standard definitions and classes of distributions. For any 
functions a{x) > 0, b{x) > 0, 

i) a{x) = o{b{x)) means lim^^oo wfl = 0; 
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ii) a{x) = 0{b(x)) means < lim^_,oo < oo; 



iii) a{x) ~ b{x) means limj-^oo = 1- 

The class Ti. of heavy-tailed distributions consists of those distributions F for which 
e*^F((ia;) = oo for all t > 0. A subclass of the heavy-tailed distributions is the 
class C of long-tailed distributions, defined by the requirement that 

x~>oo F{x) 

It is clear that, for any a G M, F{x — a) F{x). This in turn implies the existence of 
an increasing concave function h{x) < x/2 such that h{x) ^ oo as x ^ oo, see Foss 
and Zachary [T5], and 

F{x-h{x)) ^F{x). (2) 

The class S of subexponential distributions was introduced by Chistyakov in 1964 [6J. 
A distribution F on the positive half-line belongs to S if and only if 



F*2(x) ~ 2F{x). 

The class 5 is a subclass of C (see Kliippelberg [20j). By induction we also have 
F*'^{x) ~ nF{x) for any n G N. If F G 5, then for any function g{x) — >• oo as a; — >• oo, 

g{x) < X, _ 

F{x - y) 



-F{dy) — >• as a; — >• oo. (3) 

'g(x) F{x) 

This result, to be found in [llj, follows from the observation that (|3| is valid for any 
g{x) = h{x), where h{x) satisfies ([2]), and if ^ is vahd for some g{x), then it is valid 
for any gi{x) such that g{x) < gi{x) < x. For convenience we define the following 
notation. 

K^ix,r)^K{x,r) := ~ - 1, (4) 

F[x) 

r'"" F(x - v) 

Mx,r) ^ J{x,r) := \^ ,^ F{dy). (5) 

Jr F[x) 

Both these expressions have natural interpretations for, respectively, long-tailed and 
subexponential distributions. 

We note that, for F G 5, 

i) J(a;, r) — ^ F{r) as x — > oo for fixed r; 

ii) J(x, r) is monotonically decreasing in r for fixed x; 

iii) K{x^ r) ^ as X oo for fixed r; 

iv) K{x,r) is monotonically increasing in r for fixed x. 
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Lemma 2.1. Let F E C Then F E S if and only if for any c > 1 there exists a 
function h{x) = h{x, c) and a constant Xi = Xi{c) > such that, for x > Xi, 

i) h{x) is an increasing concave function ; 
a) h{x) — > CX3 as X ^ oo; 

Hi) there exists Xq > such that h{x) < x/2 for x > Xq; 
Hi) J{x,h{x)) < cF{h{x)); 
iv) K{x,h{x)) < F{h{x)). 
Proof. Let F & S and c > 1. Then there exists /c G N such that 

k = mm J G N : mf ^ > - 

\ x>j F{x) c 

Let X-i = xo = 0. 

For fixed t > k, J{x,t) F{t) < cF(t + 1). So, for r G N, we can inductively define 
an unbounded, increasing sequence: 

Xr = inf (a; : x—Xr-i > Xr-i—Xr-2 > 3, sup J{y, k+r—1) < cF{k+r), sup K{y, k+r) < F{k+r)). 
Now define 

kx/xi X < Xi, 

k + r-l + X e[xr,Xr+i),r eN. 



h{x) = 

For X G [xr,Xr+i),r G N, we have 

k + r — 1 < h{x) < k + r. 

So, _ _ 

J{x, h{x)) < J{x, k + r-1) < cF{k + r) < cF{h{x)) 

and 

K{x, h{x)) < K{x, k + r) < F{k + r) < F{h{x)). 

By construction increasing, unbounded above and concave, and, for x > Xi + 6K, 
h{x) < x/2. 

On the other hand, if F is long-tailed, and such a function h{x) exists, then by 
Proposition 2 of [4^ F is subexponential. □ 

3 Upper Bound for the Relative Accuracy 



The main goal of this section is to obtain Theorem 3.1 Our intention was to make the 



statement and proof of this theorem look very similar to Theorem 3 in [19]. However, 



the quantities defined in the statement of our Theorem 3.1 are not identical to the 



corresponding quantities in Kalashnikov and Tsistsiashvili's paper. 
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The structure of this section is as follows. First we discuss the accuracy of the asymp- 
totic approximation We construct an upper bound for P{S„ > x\v > 1). We 
introduce some notation for two common expressions associated with subexponential 
distributions, and use these to construct an upper bound for the relative accuracy 
A{x). We then introduce our test function g{x), whose existence is guaranteed by 
Lemma 2.1 and re-express the upper bound for A(x) in terms of C{x), the ratio 



between A(x) and the test function g{x). Finally we prove Theorem 3.1 



3.1 Accuracy of the Asymptotic Approximation 

It is well known (see, for example, [9]) that asymptotically 

P(^^ > x) ~ EuF{x). (6) 

In general the accuracy of the asymptotic relation ^ is poor for small values of x. 
For instance, if we consider a Pareto(5,l) distribution, and take p = 0.2, so 



F{x) 



1 for X < 1, 
x~^ for X > 1 



then, estimating the tail of the geometric sum using naive simulation, 

P{S^ > 30) = 0.00547 



E{u)P{X > 30) = 2.06 X 10"^ 

This gives us a relative error of more than 26000 in using the asymptotic expression. 
Hence the need to establish an upper bound for the relative accuracy of the asymptotic 
relation. 

Let _ 

^ ^ ^ ^ EuF{x) 
and A[a,b] = sup^<^<f, A(x). 



3.2 Upper bound for P(5'j, > x\iy > 1) 

We are now ready to construct an upper bound for P{Si, > x\v > 1). 

Choose h{x), such that F{x—h{x)) ~ F{x). We consider the total probability formula: 

P{S^ >x) = P(Xi > x)P{u = 1) + P{S^ > x\v > l)P{u > 1). (7) 

Conditional on u > 1, S^, — Xi has the same distribution as S^. Let S^, be a random 
variable, independent of {Xi} with the same distribution as S^- Denote the event 
{Si, > x,iy > 1} := B{x). We can then partition the event {Sy > x,u > 1} as follows: 

B{x) = {Xi < h{x),B{x)} U {h{x) < Xi < x - h{x),B{x)} 

U {Xi > X - h{x),B{x)}. 
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We then have 

P{S^ > x\v > 1) <P(Xi < h{x), S^> X- h{x)) 

+ P{h{x) <Xi<x- h{x), > X - Xi) + P(Xi > X - h{x)). 

We know that 

P(Su >x) = -(1 + A(x))F(x). 
P 

Hence, 

P{S^ > x\v > 1) <P(Xi < h{x))P{ST^ >x- h{x)) 

nx—hix) ^ 

+ / P{Sv >x- y)PiXi e dy) + P(Xi > a; - h{x)) 

Jh(x) 

<-(l + A(x - h{x)))F{x - h{x))F{h{x)) 
P 



/•x—h(x) 

+ -{1 + A[h{x),x - h{x)]) / F{x-y)F{dy) + F{x-h{x)). 

P Jh{x) 



3.3 An Upper Bound for A{x) 

To simphfy our expressions we recall the quantities defined in (3) and (4), and denote 
J = J{x,h{x)) and K = K{x,h{x)). We note that, since F G 5, we may choose a 
concave function h{x) oo such that both J and K converge to as x — > oo. 

Using such a function H we can now construct an upper bound for the relative accu- 
racy A(a;). 

We have 

^^-.^""^ <— ((1 + A{x - h{x))){K + l)Fihix)) + (1 + A[h{x), x - h{x)])J) 
F[x) p 

and so, 

A(a;) <(1 - p) (A(x - h{x)){K + l)F{h{x)) + A[/i(x), x - h{x)]J) 

+ (1 - p) J + (1 - p')K -{l-p){K + l)F{h{x)). (8) 

Our upper bound for A(x) should be a monotonic function Am(x), decreasing to 0, 
such that, for some constant 6, 

A(x) < A„(x) for all x>h. 

3.4 The Test Function 

As in [19] we now introduce our test function g{x). 
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We require g{x), which will depend on our choice of h{x), to have the following 
properties as x oo: 

g{x) monotonically; (9) 

max( J(a;, h{x)), K{x, h{x))) = 0{g{x)); (10) 

g{x - h{x)) r--^ g{x); (11) 



The existence of such a function is again guaranteed by Lemma 2.1 since we may 



choose g{x) := F{h{x)), with h{x) concave as in Lemma 2.1 The relation h{x 



h{x)) > h{x) — h{h{x)\ for large enough x, then implies (10). Define 

( 

C(x) = max I 0, — -— and C[a, 61 = max Cix). 

\ g{x) J a<x<b 



Now, if we set 

fi{x) = 

f2{x) = 

fsix) = 

we can rewrite (7) as 



;i - p)g{x - h{x)){K + l)F{h{x)) 

9{x) 
'l-p)g{h{x))J 

9{x) 

[1 - p)J + (1 - p^)K -{l~p){K + l)F{h{x)) 
9{x) 



C{x) < (/i(x) + f2{x))C[hix),x] + fsix). (12) 

Also, /i(x) — » 1 — p, /2(x) — > 0, and /3(x) is bounded from above. Let S = S{x) = 
suPy>a;((/i(y) + f2{y))- Hence we can find b > such that 

ifiix) + f2{x)) < 6{b) < 1 and fs{x) < ip{b) = ^ for all x>b. (13) 
3.5 The Key Result 

Following the arguments in jl9] we can now prove: 

Theorem 3.1. Let F E S, h{x) satisfy condition and g{x) satisfy conditions (8), 
(9) and (10). Then, there exists 6>0, 0<(5<1 and > 0, such that for all x > b, 

A{x) < A„(x) := Cg{x) 

where 

C = max(^-^,y^ + 6C[h{b),b]y (14) 
Proof. For a; > 6 it is clear that A(a;) < C{x)g{x) < C[b,x]g{x). 



From (12) we know that C{x) < {fi{x) + f2{x))C[h{x)^ x] + fsix). Hence, with b,5 

C[b,x] < 6C[h{b),x] + if. 



and as in (13) 
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Now, C[h{b),x] = max{C[h{b),b],C[b,x]). 
UC[h{b),x] = C[h{b),b] then 



C[b,x] < 6C[h{b),b]+^. 

l{C[h{b),x] = C[b,x] then 

C[b,x]<^^. 

This completes the proof. □ 

This result can be applied to the queueing problem described in the introduction by 
using the estimate_|T|. We have upper and lower bounds for the tail distribution 
and we call these F+ and F~ . We estimate the relative error, A^+, in using the 
asymptotic approximation P{Si, > a;) ~ E(i/)F+(a;). All the results for Theorem 1 
follow through using in place of F. We need to evaluate C{x) for h{b) < x < b, 
and we propose to do this using naive simulation. We shall bound C{x) above by 
considering a geometric sum random variable whose increments have distribution 
F+j so that 

Cix) < C^x) := (p^^^^ - l] /g{x). 



4 Applying the Result 

Some numerical estimation of the tail distribution of the geometric sum must be done 
in order to evaluate C[h{b),b]. The greater the value of 6, the more accurate the upper 
bound becomes. However, this comes at the greater computational cost of numerically 
evaluating the tail of the distribution. A compromise has to be struck between the 
tightness of the upper bound and the resources one is willing to invest in evaluating 
the tail. 

A critical part of the procedure is the choice of the test function g{x), which itself 
depends on the choice of h{x). In [19] the function g{x) was chosen as a function 
in closed form over the whole range of values of its argument. We observe, however, 
that in evaluating C[h{b),b], we know the (numerically) exact value of A(a;) in the 
range h{b) < x < b, and we wish to use this information in our choice of g{x) when 
appropriate. Our strategy in applying the result is therefore as follows. 

1) Decide what resources are available for estimating the tail distribution to a suit- 
able degree of accuracy. Given the available resources, define B as the maximum 
value for which we estimate A{B) and numerically evaluate C[h{B), B]. 

2) Determine the class of functions that will do for h{x). 

3) Estimate J{x, h{x)) and K{x, h{x)) in the range h{B) < x < B. 

4) Choose monotonically decreasing g{x), which will depend on our particular 
choice for h{x), and which may incorporate our numerical knowledge of C{x), 
such that g{x) = 0(max( J(x, h{x)), K{x, h{x)))). 
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5) If sup^>^(5(x) < 1, we take b — B, and find the corresponding value of (p. If 
suPa;>B ^{x) ^ 1 either the procedure has failed, or we must be prepared to use 
a larger value of B. 

6) Calculate C. 

Some comments on these steps will be useful. 

In step 2 when choosing h{x), there is a tension involved between the relative rates of 
decay of J{x, h{x)) and K{x, h{x)). The larger h{x) is the smaller J{x, h{x)) becomes, 
but the larger K{x, h{x)) becomes, and vice versa. We can change the rates of decay of 
J{x, h{x)) and K{x, h{x)) by scaling h{x) by some numerical factor without affecting 
the asymptotic decay rate of g{x). 

In Step 4 we will generally want to choose g{x) in order to cause the upper bound 
for the relative accuracy to decay to zero as fast as possible. We may also want to 
incorporate the information we have already calculated for C[h{B),B]. The fastest 
asymptotic decay rate for g{x) is obtained by an optimal choice of h{x). However, we 
can use the information we have gathered in calculating the numerically exact value 
of A(a;) in the range [h{h), h] by constructing a monotonically decreasing version of it, 
^Tn{x) := sup3,<2^<(,* ^{y)-i X < b*, where 6* < 6 is chosen to minimize the value of 
C. Thus, once we know the optimal asymptotic function g{x), we instead use gi{x): 

/ X f Am(a;) for x < b*, 
for some constant K chosen to make gi{x) continuous at b*. 

Hence we have a two parameters that wc can adjust, the scale factor for h{x), which 
alters the balance between the decay of J{x,h{x)) and K{x, h{x)); and the value of 
b*, which allow us, given our chosen value of -B, to minimize C, and hence tighten the 
upper bound. 

We will now show how to apply our result to Pareto and WeibuU distributions with 
various parameters. The values of J{x, h{x) were estimated using the integrate func- 
tion in R, and the naive simulation to estimate C{x) was also performed in R using 
samples of size 5 x 10^. The bounds we calculate (labelled KT bounds in the graphical 
displays to reflect the original source of this method in the work of Kalashnikov and 
Tsitsiashvili) are compared to values of the relative error that were calculated using 
a discretized Panjer algorithm with bandwidth of 0.005. 



4.1 Pareto Distribution 

We will consider Pareto Distributions of the following form: 

Fix) = [ ° for X < 1 
^ > \ 1 - for X > 1 

where a> 1. We now follow the steps above. The choice of h{x) is determined by the 
requirement that K{x,h{x)) — > 0. This occurs if and only if h{x) — o{x). We then 
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have 



X ) ~ {x- h{x)Y+^' 

The simplest form of h{x) is h{x) = x'^ for < /3 < 1. We can then choose 

g[x) = a;-i»in{a/3,l-/3)^ 

which ensures that g{x) = 0{ma:x{J{x,h{x)),K{x,h{x)))). If we want to make g{x) 
decay as fast as possible the optimal choice for h{x) will have (3 = 1/(1 + a). 

4.1.1 Example 1 

We shall take a = 2.2 and p = 0.5. We use a discretized Panjer recursion to estimate 
the tail of the distribution, and assume that our resources allow us to estimate this 
up to i? = 100 using a bandwidth of 0.005. If we follow the approach of [19] we take 
h{x) = x^/^-^ and 

We find that 5(100) = 0.786, and ^[5, 100] = 14.4. This results in C = 13.2, giving 

A(x) < 13.2x-°•^*^^ X > 100, 

as the (corrected) Kalashnikov-Tsitsiashvili bound. A small improvement can be 
made by taking h{x) = 1.70x^/^'^, with a consequent value of C = 12.9. A further 
improvement can be made, at some cost to computational time , by choosing the test 
function g{x) to be equal to A^ix), the monotonically decreasing version of the exact 
value of A(x), up to some value b* < B, and then ensuring continuity at b*. The 
best that can be obtained with our value oi B = 100 is to take h{x) = 1.14x^/^'^, and 
b* = 21.3, so that 

Amfx) for X < 21.3, 



^^^^ ^ 8.52X-0-6875 for a; > 21.3. 



This results in the upper bound 

A(x) < 8.53x-°•^^^^ X > 100. 



The logarithm (base 10) of this upper bound has been graphed in Figure 1(a) (labelled 
as Upper Bound), along with the logarithm of the numerically exact result obtained 
from the Panjer recursion (labelled Exact). 

4.1.2 Example 2 

Now consider a = 2.2 and p = 0.2. Once again we take B = 100. If we follow 
the methodology of [12] and take h{x) = x^^^'^ and g{x) = we find that 

min(n G N : 6{n) < 1) = 1085 which not only is greater than our chosen B, but 
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impracticably large in any event. No improvement can be made just by scaling h{x). 
However, if we take h(x) = 1.054a;^''^'^ and adjust g(x) to coincide with Am{x) for 
X < b* = 27.1, we arrive at 



g{x) 



Arn{x) for X < 27.1, 
^26x-o-6875 for 3, > 27.1. 



This gives the upper bound 

A{x) < 179.85x-°-6^^ X > 100, 



which is shown in Figure 1(b) 



4.1.3 Example 3 

When the Pareto distribution in question is lighter tailed the asymptotic approxima- 
tion becomes drastically less good for moderate values of the argument, and it requires 
more resources to compute numerically the tail distribution of the geometric sum for 
higher values of the argument. As our example we take a = 5,p = 0.5 and shall 
again perform this numerical exercise using a discretized Panjer algorithm, but with 
B = 50 and a bandwidth of 0.002. Proceeding as in fW], we take h{x) = x^l^ and 
g{x) = x-^/^. We find that 5(50) = 0.996, and C = 2215. This large value of C reflects 
the fact that the largest errors in the asymptotic expression for the tail of the sum 



occur very early on (as can be seen from 1(c) ). By adjusting h(x) to h{x) = lA6x^^^ 



we can reduce C to C = 1662. However, in order to avoid the initially large values of 
A(x) we take b* = 27.1, h{x) = 1.94x^/^ and 



Am{x) for X < 27.1, 
93.93X-5/6 for x > 27.1, 



which yields a dramatic improvement, giving 

A(x) < 93.7x~^/^ X > 50. 
This bound is shown in Figure l(c)[ 



4.1.4 Example 4 

For this example we return to the example of the classical Gl/Gl/1 queue as de- 
scribed in the introduction. Suppose that we have obtained an upper bound for the 
distribution of ipi as in ([T]) with distribution F+, so that P^ipi > x) < F+(x) = 
+ lx~^ for X > 1, and that we have a value of p = 0.5. We estimate an upper 
bound for the tail of the distribution of the geometric sum using a naive simulation of 
a geometric sum with increments having distribution F+ and decide we have enough 
resources to B = 80. We choose h{x) = x^^^ and g{x) = x~'^^^. We find that C = 13. 
Hence 

A{x) < 13a;-2/^ x > 80. 
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(a) Ex 1. Pareto, a = -2.2, p = 0.5, (b) Ex 2. Pareto, a = -2.2, p = 0.2, 
A„(x) = 8.53a;-°-6875 A„(a;) = 179.85x-o-6875 
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(c) Ex 3. Pareto, a = -h, p = 0.5, (d) Ex 5. WeibuU, /3 = 0.5, p = 0.5, 
A„(a;) = 93.7a;-5/<5 A„(a;) = 2.Qh2K{x, h{x)) 

Figure 1: Examples 1,2,3,5. Plot of log]^o(^(^)) logio(Au(a;)). 



4.2 WeibuU Distribution 

We will consider WeibuU Distributions of the following form: 



F{x) 



for X < 0, 
1 — e"^" for X > 



where < P < 1. Straightforward calculations show that h{x) must be chosen so that 
h{x) — o{x^'~^). Further calculations show that, if we take (log(x))^/^ = 0{h{x)), 
then J(x, h{x)) — 0{K{x, h{x)), and hence we may take g{x) — K{x, h{x)). 
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4.2.1 Example 5 



For this example we take (3 = 0.5, p = 0.5. The numerical calculations were done 
using a Panjer recursion with bandwidth 0.002, and B was taken to be i? = 100. The 
optimal function for h{x) is h{x) = (log(x))^. Applying the methodology in |19], we 
find that S{x) > 1 for x < 1660, and hence no results can be obtained (for B < 1660). 
However, by taking h{x) = 0.179(log(x))^, we obtain 

A(x) < 2.952K{x, h{x)), x > 50, 



where K{x,h{x)) = exp{^/x — y x — 0.179(log(x))^) — 1. This Weibull bound is 



shown in Figure [1(d) Because of the lack of a very sharp peak in A(x), no further 
improvement can be obtained by incorporating the numerical values of A(x) into the 
function g{x). 
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